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ABSTRACT 

We provide a consistent framework for estimating galaxy counts and variances in wide- 
field images for a range of photometric bands. The variances include both Poissonian 
noise and variations due to large scale structure. We demonstrate that our statistical 
theory is consistent with the counts in the deepest multiband surveys available. The 
statistical estimates depend on several observational parameters (e.g. seeing, signal 
to noise ratio), and include a sophisticated treatment of detection completeness. The 
JAVA calculator is freely available^ and offers the user the option to adopt our consis- 
tent framework or a different scheme. We also provide a summary table of statistical 
measures in the different bands for a range of different fields of view. 

Reliable estimation of the background counts has profound consequences in many 
areas of observational astronomy. We provide two such examples. One is from a recent 
study of the Sculptor galaxy NGC300 where stellar photometry has been used to 
demonstrate that the outer disc extends to 10 effective radii, far beyond what was 
thought possible for a normal low-luminosity spiral. We confirm this finding by a 
reanalysis of the background counts. Secondly, we determine the luminosity function 
of the galaxy cluster Abell 2734, both through spectroscopically determined cluster 
membership, and through statistical subtraction of the background galaxies using the 
calculator and offset fields. We demonstrate very good agreement, suggesting that 
expensive spectroscopic follow-up, or off-source observations, may often be bypassed 
via determination of the galaxy background with GalaxyCount. 

Key virords: methods:statistical-galaxies:general-galaxies:individual:NGC 300- 
galaxies : clusters individual: Abell 2734 



1 INTRODUCTION 

With the advent of wide-field imagers on high-quality ob- 
serving sites, there are increasing numbers of very deep, 
wide-angle studies (e.g. nearby galaxies) that are subject 
to an unknown contribution of background galaxy counts. 
Even when the same field has been observed in several pho- 
tometric bands, it can be very difficult to remove correctly 
the background galaxy counts at a deep magnitude limit. 

The issue of background contamination is likely to 
persist for many years. Even with the fine resolution (~ 
0.1") of the Hubble Space Telescope (HST) the highest 
redsh ift sources are still largely unresolved l|Quchi et al.l 
I2OO5I ). Therefore star/galaxy discrimination of such objects 
awaits the development of multi-conjugate adaptive optics 
(MCAO) capable of spatial resolution on the order of tens 
of milli-arcseconds. However, higher resolution observations 



will result in a smaller field of view (FOV), since the focal 
plane of any instrument can only be populated with a limited 
number of CCDs. A "wide-field imager" in an era of MC AO 
is un likely to exceed a few arcminutes ((Bland-HawthornI 
I2OO6I ). A smaller FOV results in poorer statistical subtrac- 
tion of background galaxies, due to larger relative uncer- 
tainty on the number counts. 

We provide a JAVA calculatoiQ for estimating these vari- 
ations in several photometric bands over a wide magnitude 
range. Our analytic approach is supported by the deep- 
est surveys to date. Our modelled magnitude limits are 
([/, B, R, I, K) (27, 29, 28, 28, 23) mag. 

We first describe our model in Section [2] and then 
present the results of the model with comparison to the 
Subaru Deep Field (SDF) in Section [3 Section g] de- 
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scribes worked examples of two practical applications of 
the model. First to a s tudy of star-counts in NGC 300 by 
iBland -Hawthorn et al] (|2005l ). and secondly to the deter- 
mination of galaxy cluster luminosity functions. Finally we 
discuss the results and present our conclusions in Section [5] 
We have assumed a cosmology of Ho = 70 km s~^ 
Mpc~^, Hm = 0.3 and JIa = 0.7 throughout, although the 
results are not very sensitive to the exact choice of cosmol- 
ogy. 



2 PREDICTING COUNTS AND VARIANCES 

The variance of galaxy number counts in a randomly placed 
cell of solid angle Q is given by, 

^nn + n^ J dnidn2Uj{9i2), (1) 

where n is the number density of galaxies, and dfii and dQ,2 
are elernents o f the solid angle il separated by an angle 9i2 
l|Peebleslll98"ol ). The angular correlation function to can be 
parametrized as, 

Lo{e) = Ae-\ (2) 

Following [Peebles! (|l975h we change variables to rewrite 
equation [T] as, 

a ^nn + {nQ,) AO^ / -j, (3) 

J jXi X2| 

where the integral on the right hand side is an integral over 
the cell scaled to unit area, and the scale factor of the cell 9c 
has been taken out of the int egral. For a sq uare cell of side Oc 
the integral equates to 2.24 (|Peeblesll 19751) . for a circular cell 
the integral equates to « 1.53 with 6c equal to the diameter 
of the cell (found by numerical integration). GalaxyCount 
also allows oblong and elliptic window functions, which are 
integrated numerically within the code. 

Thus in order to calculate the variance for any general 
observation of a given area and depth, it is necessary to know 
how the number density and the correlation function ampli- 
tude change as a function of magnitude in the appropriate 
waveband. These will now be considered in turn. 

2.1 Number density 

A large number of published galaxy number count surveys 
exist in many different wavebands. These have been com- 
piled for the UBRIK filters to yield the differential num- 
ber density over a wide range of magnitudes (18.375 ^ 
U s; 27.45, 10.45 < B < 29.35, 14.42 < 7? < 28.37, 
12.25 < / «; 28.25, 9.37 < if ^ 23.48). The list of sources 
used in this work is given in appendix I All 

Thus the number density may be computed through in- 
tegration over the required magnitude range. In practice this 
is achieved by averaging all the data in bins of width 1 mag, 
and fitting a cubic spline through the resulting averages. 
The spline is then integrated numerically. 

2.1.1 Completeness 

A very important caveat with the above method is that the 
published number density surveys have all been corrected 



for incompleteness at their faintest magnitudes in order to 
estimate the true number counts of galaxies (which is nec- 
essary to accurately derive cosmological models, etc.). How- 
ever, if the variance of galaxy counts is required to estimate 
how accurately the foreground/ background galaxies can be 
removed from an observation the above method will yield in- 
correct results. The completeness of the observations must 
be accounted for. Fewer counts will be made at the faintest 
magnitudes, and hence the contribution to the variance will 
decrease. Because the number density is very high at faint 
magnitudes, the contribution from the faintest galaxies dom- 
inates the calculation, therefore it is imperative to accurately 
account for any incompleteness. 

We define completeness as the number of galaxies de- 
tected divided by the true number of galaxies. In practice 
the completeness is usually estimated by artificially gener- 
ating galaxies of various magnitude and adding these to the 
images. The completeness can then be estimated by run- 
ning the source detection software and counting the number 
of artificial objects which have been recovered. Thus the 
completeness will depend on the source detection algorithm 
and thresholds. In particular the completeness function will 
change with the significance of the detection threshold, e.g. 
a 5(7 detection threshold (i.e. the flux of detected sources is 
> 5a above the background) will have brighter completeness 
limits than a 3a detection threshold. 

The detection threshold also determines the efi'ect of 
spurious detections, such that a higher threshold will have 
fewer spurious detections. GalaxyCount does not account 
for spurious detections, which would be few for typical de- 
tection thresholds of ^ 3a. This should be borne in my mind 
by the user if it is expected that there are many spurious 
detections in their data. 

The completeness for any observation depends on the 
telescope aperture, the combined throughput of the tele- 
scope, instrument and atmosphere, the exposure time, the 
seeing (F) and the signal to noise ratio (SNR) of the detec- 
tions, providing that the seeing is well sampled. Thus if the 
completeness function is known for any observation it may 
easily be scaled for any other observation, on the assumption 
that the shape of the curve does not change. Furthermore 
the completeness function can be modelled using a sigmoid 
function of the form, 

f ^ l_l_e«(^-ft)' 

where / — A^'moasurcd/A'truc is the completeness, b is the 
magnitude at which / — 0.5, and a (which has units of 
magnitudes"^) determines the sharpness of the step. 

It is essential to use a sigmoid completeness function 
that correctly matches the observations in question, i.e. a 
and b must be correctly chosen. This is very important be- 
cause the galaxy number counts rise rapidly at faint magni- 
tudes, and a change in the completeness function can have 
a large effect on the mean number counts and associated 
error. In Figure [1] we have compiled various completeness 
functions from the literature and our own data to show the 
spread in sigmoid slope which results from different object 
detection methods. Object detection optimised for galaxies 
generally yield a flatter slope than those optimised for point 
source detection. Similarly adaptive optics corrected images 
seem to result in a steeper slope than natural seeing images. 




-4 -2 2 4 

Magnitude 



Figure 1. Completeness functions compiled from the literature 
and our own data, all translated to have 50 per cent completeness 
at zero magnitude. The points show the data from: the Subaru 
Deep Field B band galaxy counts (open circles. [Xashikawa et al] 
\2Q04) : the author's own CMOS g band observations of a red- 
shift 0.7 cluster of galaxies (open squares); the author's own 
IRIS2 K band observations of the same field (open triangles); 
ISAAC 1.19/im narrow ba nd filter data of t he Chandra Deep 
Field South (closed circles, ICubv et al.ll2007^: NACO AO cor- 
rected observations (closed sq uares . ICrescT et al.ll2006f) : the au- 
thor' s own CMOS r band irnaging for point sources (closed trian- 
gles, |Bland^H^«thor^^er^l]|2TO^]_^^ observa- 
tions for point sources (stars. ICresci et al.ll2006l V The lines show 
the best fitting sigmoid functions and range from a = 1.4 for the 
SDF observations to a = 6.0 for the NACO AO corrected point 
sources. 



These trends are easily understood as extended sources of 
fixed magnitude will have a greater range of surface brigiit- 
ness than point sources of fixed magnitude. Tiius tiiere is a 
larger spread in the sigmoid transition for extended sources 
as some faint magnitude may have a high central surface 
brightness and thus be detected, whereas some bright mag- 
nitude sources may have an overall low surface brightness 
and thus be missed. Only sigmoid functions derived from 
extended source detections should be used with Galaxy- 
COUNT because of this difference. GalaxyCount uses a 
default slope of a — 3.0, which is typical for galaxy detec- 
tions, however we stress that to obtain accurate results the 
completeness function for the users own data should be com- 
puted. This illustrated in Figure [21 which shows the effect 
of varying the slope of the sigmoid function for B band ob- 
servations with 18 ^ -B ^ 27 mag and a fixed 50 per cent 
completeness of B—27 mag. 

GalaxyCount uses a default sigmoid with a = 3.0 
which may be used to provide an estimate the expected com- 
pleteness of any observation by scaling b appropriately. The 
b values were me asured from the SDF B, R and / complete- 
ness functions in lKashikawa et al.l (|2004l ). and converted to 
Vega magnitudes using the transformations given in sec- 
tion l2.4l ([/ was derived from the B band completeness func- 
tion), and t he SDF K band comp leteness function was mea- 
sured from iMaihara et al.l (|200ll ). In order to allow scaling 
from the SDF observations we have taken throughputs for 
Suprime-Cam BRI (actually i') from lMivazaki et al.l (|2002l ) 
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Figure 2. The effect of the sigmoid slope on the resulting number 
counts and standard deviation for 18 ^ S ^ 27 in a 5 X 5' area, 
with a fixed 50 per cent completeness of 27 mag. The solid line 
shows the number counts and the shaded area the la uncertainty 
on those counts. 



(convolved with M auna Kea atmospheric tr ansmission from 
iHook et al.ll2004l ). and for CISCO K from iMotohara et all 
(|2002l ). we have assumed the CMOS throughpu t for U. We 
also o ffer the Gemini CMOS throughputs from iHook et al.) 
(|2004 ) as a default option in the calculator (we assume that 
the u'g'r'i' throughputs are appropriate for UBRI, and we 
assume the CISCO K band throughput is appropriate for 
Gemini NIR imagers). 

Provision is made for the user to override any of our as- 
sumptions when using the calculator. Observing conditions, 
telescope aperture and instrument throughput may all be 
freely specified, or a particular sigmoid completeness func- 
tion may be input. To give an impression of the importance 
of correctly specifying the appropriate observing conditions 
etc. Figure [3] shows the effect of varying the throughput on 
the resulting 18 ^ i? 27 mag number counts and stan- 
dard deviation in a 5 x 5' area, for a 595 minute exposure, 
with SNR=3, r ~ 1", on an 8m aperture telescope. Clearly 
it is important to use the correct input for the particular 
observations of concern. 

The difference in standard deviation resulting from in- 
complete and corrected observations is shown in Figure |4j 
for a 5' X 5' FOV for a series of different exposure times and 
faint magnitude limits, with the other observing conditions 
the same as for the SDF. 

Of course the number of galaxies in the field may al- 
ready be known, and therefore GalaxyCount includes pro- 
vision for the user to enter his or her own value, as well 
as providing estimates from including and excluding incom- 
pleteness. The standard deviation a2 associated with a mea- 
sured number of galaxies, N2, counted in a particular mag- 
nitude range and area, may be scaled from the predicted 
number and standard deviation, Ni and ai, assuming the 
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Figure 3. The effect of the throughput on the resulting number 
counts and standard deviation for 18 ^ B $J 27 in a 5 X 5' area, for 
a 595 minute exposure, with SNR=3, F ^ 1", on an 8m aperture 
telescope. The solid line shows the number counts and the shaded 
area the la uncertainty on those counts. 



amplitude of the angular correlation function is the same in 
both observations, through the equation, 

^2 = ^Jpol+N^f{l~ f), (5) 

which may be combined with equation |4l to scale the results 
presented in Table [T] for different values of completeness. 



2.2 Stellar contamination 

The severity of stellar contamination will depend on the di- 
rection of the observations. Stars are not expected to con- 
tribute significantly to the num ber counts for obs ervations 
away from the Galactic plane (|Robin et al.l |2003| ) . Bright 
stars may be easily removed from the catalogues via visual 
inspectio n, radial profile fitting or automatic discrimina- 
tion fe.g. lBertin fc ArnoutsI 19961 ) . At fainter magnitudes the 
stellar contamination will not be significant. Figure [S] com- 
pares the expected contribution of star counts at a Galactic 
latitude of 6 = 90 degrees and longitudes / = 0, 90, 270 de- 
grees to the K band galaxy counts compiled from the liter- 
ature (see appendix lAll If stellar contamination is expected 
to be a problem then it can be accounted fo r using synthetic 
models such as those of I Robin et all l|2003lfl 



2.3 Amplitude of the angular correlation function 

The amplitude of the angular correlation function decreases 
at fainter magnitudes, since the average redshift of the sam- 
ple is larger and hence the volume being surveyed is larger, 
and thus the distribution of galaxies is approaching homo- 
geneity. 
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Figure 4. The expected number, standard deviation, and stan- 
dard deviation divided by square root of the number (i.e. the 
non-Poissonian component of the standard deviation) of galax- 
ies in surveys of different magnitude limits, with a 5' X 5' FOV. 
The thick black line is for the total number of galaxies, the dot- 
ted and dashed lines include estimates for the incompleteness of 
the observations for exposure times of 595, 300 and 10 minutes, 
with SNR=3, r 1", with an 8m aperture telescope with 36% 
throughput and deviate from the thick black line in that order. 



The functional dependence of A on magnitude 
has been reported to be monotonically dechning (e.g. 
iMcCracken ct al.' 2001; Wilson 20031: ICoil et al."2004 !) or to 
flatten off at the faintest ma gnitudes (e.g. Brainerd fc Smaill 



1 19981 : IPostman et al.l Il998h . Both possibilities have been 
modelled in GalaxyCount, and the user may choose be- 
tween them. The first case was simply modelled by a least 
squares fit to \ogA vs. magnitude for data compiled from the 
literature. However, to allow the models to be extrapolated 
to very faint magnitudes the second case has also been mod- 
elled, which allows for a fiattening of A at faint magnitudes. 
The sources for the measurements of the angular correlation 
function used in this paper are listed in appendix I A2I 

The amplitude of the angular correlation function has 
been modelled using theoretical ar guments describing the 
evol ution of clustering, in particular lEfstathiou et ahl (Il991h 
and iPeeblesI ([1980). The method takes a parametrized evo- 
lution of the two-point spatial correlation function and 
converts it to the expected dependence of A on magni- 
tude. The two-point spatial correlation function is usually 
parametrized as, 
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Figure 5. The modelled star counts from iRobin et al. I ll2003l 'l 
compared to galaxy counts compiled from the literature shown 
by the points. 



(6) 



The scalinR of ^(r) with redshift, z, may be parametrized as 
(jGroth fc Pee bles 1977,), 



where 

h{z)^{l + z)-<-''+^\ 



(7) 



(8) 



where r is the proper distance and e describes the evolution 
of the clustering, with 

{0 - clustering fixed in proper co-ordinates 
— 1.2 - clustering fixed in comoving co-ordinates (9) 
0.8 " prediction of linear theory. 

In all cases with e > —1.2 the comoving galaxy correla- 
tion length decreases monotonically with z. GalaxyCount 
uses e = —1.2, 0, 0.8, 3 as selected by the user; e = 3 was in- 
cluded to model the strong evolution of clustering observed 
at faint magnitudes (ICoil et al.ll2004l ). 

For small angles u!{d) may be related to ^(r) as 
lEfstathiou et al.lll99ll ). 



^ r[(7-l)/2] A' 

"W-^ r(7/2) 

where F is the complete gamma function and A' is, 



A' 



'dz 



and, 
5W = 



h{z) 



dl~\z)idr{z)/dzy 



(10) 



(11) 



(12) 



Thus, A may be computed from equation [TO] by assum- 
ing a particular cosmology and adopting a form for dN/dz. 
We have adopted an empiri cal relation for dN / dz as a func- 
tion of magnitude given bv lCoil et all l|2004 ). 



d/V 
d7 



2/2 



(13) 



where zo = —0.84 -1- 0.05/. Other bands were transformed 
to / usin g the z = colours f r om th e synthetic stellar pop- 
ulation of iBruzual fc Chariot! (120031 ) . assuming an instanta- 
neous burst of star-formation at z = 6.4 and solar metallic- 
ity. 

Note that the empirical form used for dN/dz as a func- 
tion of magnitude is not applicable for bright magnitudes. 
Therefore we use a linear relation for magnitudes brighter 
than the point at which the slope of the theoretical curve 
and the best fitting linear relation are equal. Although aes- 
thetically more pleasing this in fact makes little difference to 
the results which are dominated by the much higher num- 
ber counts at faint magnitudes. The models are compared 
to observations (listed in appendix I A2|) in Figure |S] In most 
cases the monotonically decreasing model seems to be the 
best fit to the data. 

A subtle issue is that the fits to as a function of mag- 
nitude are calculated for a median magnitude. In the sources 
that use limiting magnitudes we have converted to median 
magnitudes using our compiled number counts within the 
magnitude limits. However, A is required over a range of 
magnitudes for use in equation [T] Therefore to calculate 
between magnitude limits we substitute n per magnitude for 
n in equation [1] and then integrate over magnitude between 
the limits in question, using the appropriate value for A at 
each point. 

The effect of the chosen model of A is shown in Figure[7] 
for a linear fit and a model with e = 3. The calculations were 
performed in the / band, with a fixed bright magnitude of 
/ = 18, a FOV of 25 sq. arcmin., and a fixed exposure time 
of 300 minutes. The effect of the flattening of A at faint 
magnitudes is to increase a. 



2.4 Magnitude systems 

GalaxyCount uses data compiled from a large number of 
heterogeneous sources, listed in appendix lAll to calculate 
the number density over a large range of apparent magni- 
tudes. These different sources comprise various filter and 
magnitude systems, and have been tran sformed to a com- 
mon system of 'Kron' total magnitudes (jKron| |l98d ) in the 
Ve ga system (B is John son and R and I are Kron- Cousins, 



Metcalfe et"alll200lf ). This system was chosen as most 



where dA is the angular diameter distance. 



of the data were already in such a system, or a similar sys- 
tem, and therefore required the least amount of transfor- 
mation. However, we recognise that this system may not 
be the choice for many future astronomical systems, given 
the widespread use today of SPSS u'g'r'i'z' filters, Pet- 
rosian magnitudes (PetrosianI Il976h and AB magnitudes 
l|Fukugita et al.ll 1999 )7 Thus in order to increase the capabil- 
ity of GalaxyCount we have incorporated transformations 
to the u'g'r'i'z' system, and list below the transformations 
to this and the AB systems. These transformations are only 
approximate, since galaxies display a large range of colours, 
surface brightness profiles, spectra, etc. and thus no single 
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Figure 6. The amplitude of the angular correlation function as a function of median apparent magnitude, compiled from the sources 
listed in appendix IA2I The models described in section [2.31 are shown. In most cases a linear model fits the data very well. 
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Figure 7. The effect of the choice of model for A. The soUd Une 
shows a model for a Unear fit to logA vs. magnitude; the dotted 
line shows a model with e = 3, as described in the text. The effect 
of the flattening of A at faint magnitudes is to increase cr. 

transformation is appropriate for all galaxies. Of course the 
user may simply compute the transformations appropriate 
to his or her own data, and input the transformed magni- 
tudes. 

Conversions to SDSS u'g'r'i'z' magnitudes were calcu- 
lated using the z = colours derived from a model spectrum 
of a passively evolving, solar mctallicity galaxy, formed at 
z = 6.4, from the libraries of iBruzual fc Charloti l|200l ). The 
z = colours are assumed for all galaxies, and are listed be- 
low, 

u (7 4-0.1 
g' 0.3 

r'Riii-f-0.13 (14) 
i' ^ J + 0.1 
z' ^1-0.2. 

For comparison to the SDF (see the next section) the 
following conversions to AB magnitudes were used, 

Uab~U+ 1.0 
-Bab ~ -B- 0.1 

i?AB~-R + 0.2 (15) 
7ab ~ / + 0.5 
Kab~K + 1.9. 

For comparison to the Hubble Deep Fie lds (sectioH31l 
we us e the following conversions given by iMetcalfe et al.l 

B F450 + 0.1 
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Figure 8. The effect of the window function on the standard 
deviation for a constant area of 25 square arcminutes. 



i?«F606 - 0.1 (16) 
/ ^ F814. 



2.5 The effect of the window function 

GalaxyCount allows a choice of window functions so dif- 
ferent detector shapes and illumination patterns may be 
modelled. The choices are square, circular, rectangular and 
elliptic. Square and circular windows are treated separately 
from rectangular and elliptic window functions for the pur- 
pose of increasing the speed of the calculations. In Figure [8] 
we show the effect of the window function on the resulting 
galaxy counts and standard deviations for observations in 
the B band with a field of view of 25 square arcminutes in 
all cases. The differences can be understood in terms of the 
different allowable separations between galaxies accommo- 
dated by the windows. 



3 RESULTS 

In this section we describe the general characteristics of the 
model via a comparison to data from the SDF. We demon- 
strate that the model reliably predicts accurate standard 
deviations for the associated number counts. However the 
number counts are not always reliably predicted by the 
model, particularly in the B band. This is shown to be a con- 
sequence of the large discrepancy between published num- 
ber densities at faint B band magnitudes. We also present 
general results of the model for different magnitude limits, 
filters, and areas. 

Note well that the model works in the Vega magnitude 
system. Therefore we have converted to the SDF AB magni- 
tudes system using the transformations given in equation llSI 
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3.1 Characteristics of the model and comparison 
to the Subaru Deep Fields 

Number counts and standard deviations have been com- 
puted from the SDF by random selection of square fields 
wholly contained within the SDF (and avoiding regions near 
bright stars and bad columns) and counting the number 
of galaxies. Following iKashikawa et alj (I2OO4I ) we distin- 
guish stars from galaxies using the following conditions: if 
mag ^ 20 and starclass< 0.99, or if 20 > mag ^ 24 and 
ln(isoarea x pix) ^ (24 — mag) /150, or if mag > 24, then 
we consider the object to be a galaxy, where starclass is the 
SExTRACTOR star/galaxy classifier, isoarea is the isophotal 
area above the SExtractor analysis threshold and pix is 
the pixel scale. For the discussion below 100 2 x 2' square 
fields were randomly chosen in each band with a fixed bright 
magnitude limit of 20 mag in each band. Integral counts are 
considered in all cases. 

GalaxyCount can be used to estimate the number 
counts of objects within the FOV and magnitude range, from 
which the standard deviation will be derived, or if the num- 
ber counts of galaxies are already known they can be input 
directly into the calculation. We present the results of three 
models with different methods of arriving at the number 
counts. All models use e = 3 for the variation of the am- 
plitude of the correlation function which produces a higher 
standard deviation at faint magnitudes than a log{A) oc m 
model, because A becomes fiatter at faint magnitudes for 
the e models (see equation [T] and Figure [T}. 

The first model uses number counts taken directly from 
the SDF. This allows us to examine the predicted standard 
deviation in isolation from the predicted number counts. The 
results are shown in the left hand panels of Figure [9] The 
top panel shows the R band SDF counts and the model. 
The grey area shows the la uncertainty around the number 
counts, and the hatched area shows the 2a uncertainty. The 
open points are the SDF mean number counts and the error 
bars are la. The middle panel shows the variation in the 
total uncertainty, a with limiting magnitude. The bottom 
panel shows the variation in the non-Poissonian component 
of the uncertainty. Whilst the model provides a reasonable 
match to the observed data for the total uncertainty, the 
non-Poissonian uncertainty is underestimated. The reason 
for the disagreement is that although n has been exactly 
specified, the variation of n with magnitude has not and 
therefore an average value for A must be used leading to the 
inaccuracy. 

The second model has number counts calculated from 
published galaxy number density surveys. The results are 
shown in the middle panels of Figure |9l and the same gen- 
eral trends are seen as for the first model. At the faintest 
magnitudes this model over predict the expected number of 
galaxies. This is to be expected since the published number 
densities have been corrected for incompleteness, whereas 
our SDF measurements have not. The over-prediction of the 
number counts contributes to the standard deviation being 
over-estimated. 

In the third model, shown in the right hand panels 
of Figure |9] the incompleteness of the SDF observations 
has been taken into account to avoid the over-prediction of 
galaxy counts at faint magnitudes, using a sigmoid with a 
slope of a = 1.4, which was found by fitting to the complete- 



ness functions from IKashikawa et all (|2004l '). This model 
matches the data extremely well, both in terms of the total 
uncertianty and the non-Poissonian uncertainty. 

The models were also tested for in the B and i' bands, 
and show very similar results. However, for the number 
counts at faint limiting B and i' band magnitudes the model 
overcompensates for the incompleteness despite the fact that 
the incompleteness fun ction is ta k en dir ectly from the SDF 
observations of Kashikawa et al] (|2004l ). Consequently the 
predicted uncertainties also fall below the measured values. 

These discrepancies may be due in part to uncertainties 
in correcting for both the magnitude systems and the differ- 
ent filters used by the SDF and the collated number counts. 
Furthermore, the number counts are calculated from a cubic 
spline fit to an average of compiled number density surveys 
from the literature. The B band SDF counts are systemat- 
ically higher than the other surveys which reach faint mag- 
nitudes, and thus the average counts are already lower than 
the SDF counts before the incompleteness correction is ap- 
plied. However, this still leaves the problem that the uncer- 
tainties in number counts are under-predicted by the model. 
We have shown in Figure [5] that the model will predict the 
uncertainties very accurately if the number counts as a func- 
tion of magnitude are correct. Therefore it seems that the 
B band SDF counts differ from the average counts by more 
than would be expected from cosmic variance. This is not 
to say the B band SDF counts are wrong, but rather that 
there are systematic differences between published number 
count surveys that are larger than would be expected from 
cosmic variance, most likely resulting from uncertainty in 
correcting for incompleteness and subsequent extrapolation 
to the standard one square degree for which number densi- 
ties are quoted. This is illustrated in Figure fTUl which shows 
number density vs. magnitude. It is clear that the scatter 
dramatically increases at B ^ 25, which is unlikely to be a 
real effect due to cosmic variance. In comparison the scatter 
in the R band counts is much more uniform at faint mag- 
nitudes. The expected uncertainty due to cosmic variance 
is also shown in these plots. The grey curves show the ex- 
pected cosmic variance from a survey of 1 deg'^; the hatched 
curves show the expected cosmic variance from a survey of 
1.17 X 10"^ deg2 (the area of the Hubble Deep Field) extrap- 
olated to 1 deg^. Both simulations are assumed to be 100% 
complete, an additional small correction would be needed 
for incomplete surveys. The published number densities dif- 
fer within cosmic variance at bright magnitudes but have a 
larger scatter at faint magnitudes. Whilst this is markedly 
more pronounced in the B band it is also true in the R band. 
The extra scatter at faint magnitudes can be attributed to 
systematic errors in accounting for the incompleteness of the 
surveys, or to systematic errors in consolidating the differ- 
ent magnitude systems of each survey. Note that an accurate 
number density is crucial to correctly calculate the associ- 
ated uncertainties, and thus the average presented here may 
fail to achieve this if systematic errors really are present. 
We note though that the user is free to specify the number 
density exactly, and so may choose to rely on a particular 
survey of his or her own choice e.g. the Hubble Deep Fields, 
or to input a number density from his or her own survey. 
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Figure 9. Comparison of tlie model to the SDF R band. The model number counts have been calculated for an e = 3 model. The left 
hand panels show the data matched exactly to the SDF number counts. The middle panels show a model with number counts calculated 
from number density surveys (see appendix lAll l. The right hand panels show a model with number counts calculated from published 
surveys, but corrected for the incompleteness of the SDF observations. The top panels show the SDF mean number counts and the 
associated stand deviation, as open points with error bars. The grey area represents the predicted la uncertainty from the model, and 
the hatched area the 2a uncertainty. The middle panels show the increase in c as a function of limiting magnitude, and the bottom 
panels show the variation in a/\/T[, i.e. the non-Poissonian component of the variance. 



3.2 General results 

Table [1] shows the results for a selection of magnitude ranges 
and areas. These values were computed using the e = 3 
model, and assume 100% completeness at all magnitudes. 
The variances may be easily scaled to any number counts 
within the same magnitude range and area, assuming that 
the amplitude of the correlation function does not change, 
equation [5] After appropriate scaling, equation |4] may be 
used to estimate the completeness if the number of galaxies 
is unknown, e.g. if star-galaxy discrimination is not possible. 

We have explicitly demonstrated the reliability of the B, 
R and I band predictions above, via comparison with the 
SDF. However, the U and K bands have not been checked 
since there are no sufficiently large and deep surveys to per- 
form a satisfactory random counts-in-cells test. The K band 



has good statistics on the number counts and the angular- 
correlation function, so there is no reason to suspect that the 
results will be very wrong but until they are checked against 
observations we cannot know for certain. The U band is still 
more uncertain since the variation of A with magnitude is 
unknown, and so the B band function has been assumed, 
and normalised to the U band. 

4 WORKED EXAMPLES 
4.1 NGC300 

We now provide a worked example to emphasise the im- 
portance of accurately determining the background galaxy 
count s in a given photometric band. iBland-Hawthorn et al.l 
(|2005l ) obtained deep observations of the stellar population 
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Figure 10. The number density of galaxies vs. magnitude, with B on the left and R on the right. The grey shaded curves show the 
number density and la uncertainty expected from cosmic variance. The hatched curves show the cosmic variance as measured in a survey 
with the the area of the Hubble Deep Field, and subsequently extrapolated to 1 deg^ . 



in the outer disc of NGC 300, a late-type spiral in the Sculp- 
tor group. They used the Gemini Multi-object Spectrograph 
(GMOS) on the Gemini South 8m telescope in exceptional 
conditions (0.6" FWHM seeing). At a Sa point source detec- 
tion threshold of r' = 27.0 mag, they were able to trace the 
stellar disc out to a radius of 24', or 2.2 x R25 where R25 is the 
25 mag arcsec"^ isophotal radius. This corresponds to about 
10 scale lengths in this low-luminosity spiral {Mb ~ —18.6), 
or about 14.4 kpc at a distance of 2.0 Mpc. 

These autho rs, and others (e.g. iFerguson et al.l l2005l : 
llrwin et al.ll2005l ). demonstrate the profound importance of 
usi ng star counts for trac i ng th e outer reaches of galax- 
ies. iBland-Hawthorn et al.l (|2005| ) reach an effective surface 
brightness of 30.5 mag arcsec"^ {2a) at 55% point source 
completeness which doubles the known radial extent of the 
optical disc. These levels are exceedingly faint in the sense 
that the equivalent surface brightness in B or V is about 32 
mag arcsec"^. 

In the NGC 300 study, the authors found no evidence for 
truncation of the stellar disc such that the disc may extend 
even further. But this conclusion is sensitive to the assumed 
background galaxy counts. The disc can be forced to truncate 
at smaller radius if the background galaxy counts are close 
to or exceed 200 gals arcmin"^. However, the authors found 
in favour of a lower background (130 gals arcmin"'^) by as- 
suming that the outermost field is dominated by background 
galaxies. 

Note that in order to statistically subtract background 
galaxies it is often necessary to observe "blank-sky" fields, 
for an equal amount of time as the field of interest. The 
JAVA calculator obviates this expensive requirement (which 
in practice is not often met). Provided that the complete- 
ness of the target observations can be calculated, the deep 
galaxy surveys incorporated into the calculator can be used 
to subtract the contaminating galaxies and to estimate the 
associated error. 

We now use GalaxyCount to predict the expected 
galaxy counts in r' for a 5.5' x 5.5' field. This takes into ac- 
count the seeing (0.6"), SNR=3cr, the telescope aperture 
(8m), total system throughput=39%, and the exposure time 



(8100s). For our values (in this case we use R = r' — 0.2 < 
26.8, which is appro priate for stellar populations i n the outer 
discs of galaxies, see lBland-Hawthorn et al]|2005h . we deter- 
mine 85 ±11 galaxies arcmin"'^, or 2560 ± 110 galaxies over 
the GMOS FOV, where the errors are la, calculated using 
an e = 3 model. This is equivalent to a total completeness 
of ~ 56% for R 5C 26.8 for galaxies as predicted by the cal- 
culator. 

We now consider the impact of the background galaxy 
counts on the luminosity p rofile derived from the stellar 
counts. The primary aim of iBland-Hawthorn et al] (|2005l ) 
was to obtain deep observations of the outer disc of 
NGC 300, and hence only a section of the disc was stud- 
ied, located approximately 9kpc SW of the galaxy nucleus 
on the semi-major axis. Star counts were made in arcs of 
fixed width centred on the galaxy nucleus, thus the area of 
each segment is different. Consequently the uncertainty on 
the star counts and on the background galaxy counts will 
change as a function of radius. 

Figure [11] reproduces the luminosity profile of NGC 300 
star counts. There is no evidence for a truncation of the stel- 
lar disc; the star counts decline exponentially to the limits 
of the data, which is eq ual to about 10 scale lengths. This 
confirms the findings of iBland-Hawthorn et aL* (2005), but 
with the added assurance that the errors introduced through 
the subtraction of an estimated or a theoretical number of 
backg round galaxies do not signifi cantly affect the conclu- 
sions. [Bland^Hwthor^^et^^Ll ([2003) found they could gener- 
ate a truncation of the disc if the background contains ~ 220 
galaxies arcmin"'^, but this is ruled out at high significance. 

4.2 Galaxy cluster luminosity functions 

Another application of GalaxyCount is the statistical re- 
moval of foreground and background galaxies from clusters 
in order to compute the properties of the cluster members. 
Here we provide an example of such use in the computa- 
tion of the bj band luminosity function (LF) of cluster Abell 
2734. This cluster was selected from the two degree field 
galaxy redshift survey (2dFGRS) study of the cluster LF 
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Table 1. Galaxy counts and standard deviations for a range of filters, magnitudes and areas. These values were computed using the 
e = 3 model and assume 100% completeness at all magnitudes. The standard deviations may be scaled for incompleteness or higher 
number counts using equation [5] 

Counts and standard deviation within the given area 
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(|De Propris et al.ll2003l ). based o n its hish membe r ship o f 
127 galaxies with the selection of lOe Propris et al.l l|2003t ). 
The cluster has a redshift of cz = 18646 km s~^ and velocity 
dispersion of cr = 1038 km s~^. 

We have calculated the LF in three ways. First we se- 
lect all galaxies within a box of 4Mpc centred on the cluster 
having redshifts within ±3(t of the velocity dispersion. This 
yields an accurate catalogue of cluster members from which 
the cluster LF is computed. Secondly we count all galaxies 
within the 4Mpc box regardless of their redshifts, as if the 
redshifts were unknown. We then use GalaxyCount to es- 
timate the number of interlopers as a function of magnitude 
and compute the LF based on these corrections. Finally we 
select offset fields around the cluster from which to measure 
the background galaxy counts to compare with the stan- 
dard practice of statistical subtraction of foreground and 
background galaxies. 

In order to use GalaxyCount to model the 2dFGRS 
counts the completeness limit must be accurately mod- 
elled. The 2dFGRS has a magnitude hmit of bj « 19.45 
mag, but varies slightly as a function of position. There is 
also a redshift complete ness that is a function of magni- 
tude. [Colless|^et|al] (|200ll ) provide a model for the complete- 
ness function, and software is available from the 2dFGRS 



web page^ to compute the necessary parameters to deter- 
mine the final completeness as a function of magnitude. We 
match our sigmoid model to their completeness function, as 
shown in Figure [121 Finally we convert B — bj -1-0.2. It is 
worth re-emphasising that accurate completeness functions 
are crucial for the success of modelling expected counts with 
GalaxyCount since the number counts at faint magnitudes 
are so dominant. Note though that accurate completeness 
functions can be readily computed from the imaging data in 
question. 

The offset fields were chosen at 1 and 2 degree displace- 
ments in RA and dec around cluster centre giving 24 fields in 
total. Computing the LFs for a large number of offset fields 
allows us to determine the relative merits of using Galaxy- 
Count against using offset fields, since different offset fields 
can give very different results. Note that there are two meth- 
ods in which the offset LFs may be computed: by using the 
total raw counts in the cluster and offset fields, or by count- 
ing only those galaxies for which redshifts are known. The 
former method mimics how the technique would usually be 
employed in standard observational practice, since statisti- 
cal subtraction is only performed in the absence of spectro- 
scopic redshifts, hence these are the results we report. The 
latter method accounts for the redshift incompleteness of the 

^ http:/ /www. mso.anu.edu.au/2dFGRS/ 
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Figure 11. The luminosity profile of NGC300. The top panel 
shows the raw counts as the open symbols, with the associated 
error (calculated as the square root of the counts), the solid line 
represents the expected contribution to the raw counts from back- 
ground galaxies, and the shaded region represents the uncertainty 
on this number as predicted by GalaxyCount. In the bottom 
panel the closed symbols represent the corrected star counts, and 
the combined errors. 



spectroscopic LF. Since incompleteness has been accounted 
for in the GalaxyCount derived model we also tested the 
latter technique. The results were not significantly differ- 
ent since the magnitude limit of the 2dF input catalogue is 
bright enough that observations are not much affected by 
the magnitude dependent redshift incompleteness. 

The LPs were com p uted follow i ng th e method described 
in lEllis fc Jonej (|2004l '). ISchechted l|l976h functions were fit 
for all galaxies bj ^ 19.4384, i.e. the magnitude li mit of 
the 2d FGRS at this position, through minimisation of ICashI 
statistic. The resulting fits to the characteristic mag- 
nitude, M* , and the faint end slope, a, are shown in Fig- 
ure [T31 The GalaxyCount derived Schechter function is 
inside the 68 per cent confidence limit of the best fitting 
spectroscopic results. Only one offset field is far inside the 
68 per cent confidence ellipse, with another 5 offset fields just 
inside, and another 18 further out. Thus we expect Galaxy- 
Count to be more reliable than an offset field ~75 per cent 
of the time. Forty per cent of the offset fields lie outside 
the 90 per cent confidence limits. Note that in general it 
is very difficult to get time allocation committees to grant 
more than a single offset field per cluster, and there is no 
way of knowing a posteriori whether a 'good' offset field has 
been observed. 

Another way of comparing the results is calculate the 
reduced difference between the statistically subtracted 
LFs and the spectroscopic LF. These results are presented 
in Figure 1141 which shows a histogram of the for the 
24 offset fields, with the value for GalaxyCount overlaid. 
GalaxyCount is placed in the upper 25 per cent of the 
offset fields, and occurs near the peak of the distribution. 



Figure 12. The 2dFGRS completeness as a function of magni- 
tude at the position of A2734 is shown by the solid line. The 
dotted line is our sigmoid representation. 
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Figure 13. The best fitting Schechter functions as described in 
the text. The contours show the 68 and 90 per cent confidence 
limits of the spectroscopic LF. 



i.e. it accurately refiects the expected average result that 
would be obtained from subtracting a large number of offset 
fields. 

Thus GalaxyCount provides an accurate and efficient 
way to estimate the contribution of foreground and back- 
ground galaxies to cluster LFs. In many cases expensive 
spectroscopic follow-up or off-source background observa- 
tions may be by-passed and replaced with background esti- 
mates from GalaxyCount combined with an accurate de- 
termination of the completeness of the observations. 
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Figure 14. Histogram of tlie difference between the statis- 
tically subtracted LFs and the spectroscopic LF. The value for 
GalaxyCount is overlaid and is in the upper 25 per cent of the 
offset fields. 

4.3 Comparison to the Hubble deep fields 

The Hubble Deep Field North (HDF-N) has higher galaxy 
n umber counts than th e Hubble Deep Field South (HDF- 
S, iMetcalfe et al]|200ll ). Figure [15] compares the difference 
between the two fields to the expected standard deviation 
as predicted by GalaxyCount, with the number counts 
matched to the average of HDF-N and S. The two fields 
are within la of each other when using the e = 3 model. 
Thus it seems that cosmic variance is enough to account for 
the difference between the two fields. GalaxyCount will 
be updated to include counts from the Hubble Ultra Deep 
Field when these become available. 



5 DISCUSSION 
5.1 General results 

We have presented a model to predict the galaxy counts, 
variances and completeness of observations as a function 
of magnitude, waveband and area. The models have been 
compared to observed number counts and densities calcu- 
lated using the SDF. It was found that the model will accu- 
rately predict the standard deviation on a given number of 
galaxies. Hence for any application for which the number of 
galaxies is known a priori the calculator should accurately 
determine the appropriate cosmic variance. 

For observations (such as star counts) for which the 
number of galaxies is not known, more caution must be ex- 
ercised in using the models. The derived standard deviation 
is strongly dependent on the number of galaxies. However at 
faint magnitudes an accurate estimate of the number counts 
can be difficult to determine since the scatter in values pub- 
lished in the literature varies more than the predicted cosmic 
variance (for example in the B band). The systematic errors 
introduced into the number counts are propagated into the 
resulting standard deviation. In practice this means that the 
statistical subtraction of galaxy counts may be dominated 
by systematic uncertainties rather than cosmic variance in 
the B band. However, for brighter magnitudes where com- 
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Figure 15. The difference between the Hubble deep fields N 
(squares) and S (circles), compared to the expected standard de- 
viation (l(T=grey shaded area, 2o'=hatched area). 

pleteness is not an issue the predicted number counts are 
reliable. 

The problem may be alleviated with future, more pre- 
cise determinations of the faint end number counts. Until 
such data is available it is up to the user to choose wisely 
the number counts which he or she wishes to use and to be 
aware of the possible systematic uncertainties introduced at 
faint magnitudes. 

The e = 3 model was found to predict the standard 
deviation accurately. This model has a monotonic decrease 
in log A with magnitude at bright magnitudes, but flat- 
tens off at faint magnitudes. Such a tre nd has previously 
been observed in several studies (e.g. iBrainerd fc Small! 
1 19981 : IPostman" et al.l ll998h. although many other studies 
find no such tren d (e.gjMcCrackcn ct al. 200^ : IWilson|[2003l : 
ICoil et al.|[2004h . Figure |6] shows that for the data compiled 
from the literature a simple monotonic model seems to be a 
better fit. However, we are extrapolating A to fainter mag- 
nitudes than the measurements, thus it may be that some 
fiattening is required at very faint magnitudes, beyond the 
limit of the data. Note that the counts at faint magnitudes 
dominate the statistics, thus we are testing the currently un- 
known functional form of the relation at these magnitudes. 

An e = 3 model implies a strong evolution in the clus- 
tering of galaxies since 2 = 1. However, we caution that the 
physical interpretation of the model should not be taken too 
literally, since the models do not fit A at all magnitudes, 
rather the model was used allow a flattening of A at the 
faint end. Other studies have reported that a high value of 
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Figure 16. The factor by which the total standard deviation 
dominates over the Poissonian error, a/\/N, as a function of R 
band limiting magnitude for various field sizes (labelled in the 
plot in units of square arcminutes). 

e is required to fit the faint end of the A scaling relation 
jMcCracken et alJlioOObl : iMcCracken et al ] |200ll : ICoil et alJ 
|2004| ). but note that these models cannot fit both the bright 
and the faint ends of the relation simultaneously. A precise 
determination of the value of A at faint magnitudes awaits 
future measurements. We remind the reader, however, that 
the models used in GalaxyCount can accurately repro- 
duce the measured standard deviation of galaxy counts in 
the SDF, despite the uncertainty in the precise functional 
form of A. 

GalaxyCount has a demonstrated practical use as 
shown by our reanalyses of the star counts in NGC 300, 
and the LF of Abell 2734. The calculator will prove use- 
ful for any observations subject to an unknown background 
or foreground of galaxies, and in many cases should obviate 
the need for expensive 'blank-sky' observations. 

5.2 When is cosmic variance significant? 

An important consideration is when does cosmic variance 
need to be accounted for? In some circumstances the Pois- 
sonian component of the variance will dominate over cosmic 
variance due to clustering, whereas in others the excess un- 
certainty will be very important. This has been addressed 
throughout the paper by presenting a/\^, the factor by 
which the total standard deviation is greater than the Pois- 
sonian standard deviation. This factor increases with the 
area and depth of the survey (see equation [l}. We show the 
form of the dependence in Figure [16] for the R band predic- 
tions as calculated in Table [T] 

The non-Possonian contribution quickly becomes im- 
portant as the surveyed area or the depth of the survey in- 
creases; the exact point at which it becomes non-negligible 
depends on the combination of these variables and the ac- 
curacy needed. For example for a survey reaching a limit- 



ing magnitude of i? 26 mag, Poissonian variance will be 
within a factor of 2 of the true variance only if the area is 
< 10 arcmin^. Other conditions can easily be tested using 
the values given in Table [1] 

5.3 Future work 

GalaxyCount will be kept up to date with the inclusion 
of the latest number densities and angular correlation func- 
tions as they are published, including the Hubble Ultra Deep 
Field number counts. In addition it will be extended to in- 
clude predictions for specific targets such as Ly-Q sources at 
high z, specific observing conditions such as HST/ JWST 
observations with a fixed point spread function, and ground- 
based AO-corrected observations, etc. An extension is also 
planned to include predictions of star-counts as a function 
of co-ordinates, magnitude and FOV. All future upgrades 
will be available from the web-site. 
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APPENDIX A: SOURCES OF DATA FOR 
NUMBER COUNTS AND ANGULAR 
CORRELATION FUNCTIONS 

Al Number density 

The expected number counts of galaxies in the UBRIK 
bands has been calculated using a large number of pub- 
lished surveys of galaxy number density, which have 
been comp ile d from the follow i ng pa p ers: U band, 
iKod lll986l). iGuhathakurta et all (Il990l1. iMaddox et all 
' l990l ). ISoneaila et all (Il990l). [Jones et all lll99C 



lArnouts et al. I l|200ll). iMetcalfe et al.l 1120011). ICapak et al.l 

f; B ban d. iKronI ll 19781) , IJarvis fc Tvson (198Jl 
(Il986l). IPetersoii et alT ll 19861) , Stevenson et'aLl 
, iTVsonI lll988l) iHevdon-Dumbleton et"^ il98^ 
iJones et all (Il99ll). iLillv et all lll99ll). ^Metc alfe et all 
lll99ll). iMetcalfe et all (Il995l). iBertin fc D cnncfcld lll997^^ 
IMetcalfe et al.' l'l998'). Arnouts et al' 1*1999). Arnouts ot_al] 
ll200ll). [Huang et al. (2001.), Kiimmel fc Wagner (2001 



(l200ll ). 



IMetcalfe et all (I2OOII). lYasuda et al.l .. _ _ . . . . 
( 2003f ). ICapair et al.l (|2004h . as well as the Revised Shap- 
ley Ames and the Zwicky cat. both from the D urham 
galajcy counts w eb pages; R band. Couch fc Newelll 



Hall fc Mackavl (IT984). Jnfante et al.1 (Il986l). IKocI 
m et al.1 (ll986D. lYee fc GreenI (Il987l) 
Jones et al.l (ll99ll).lMetcalfe et al.1 (Il99ll) iPicard 
I), ICouch et al.lTl99a). ISteidel fc Hamiltor] (1 19931) 
iDriver et all (ll994D. IMetcalfe et all (1 19951), ISmail et all 




|ll995l). iBertin fc Dennefeldl (Il997l) iHogg et al.1 (Il997l) 



[Metcalfe et all (Il998l) . lArnouts et "ah! (1 19991) . lArnouts etal] 
(2001), Huang et al.' ('2001'), 'Kiimmel fc Wa gneil (|200ll) , 
Metcalfe et al. ( 2001) Capak et al, (2004i); / band. 



Hall fc Mackavl (Il984l), IKooI ((19861) iTvsonI (Il988l) 



Lillv et al.1 (ll99lD. IPriVer et al l (Il994h 



Casertano et al.l 



.19951). iDriver et al.l ( Il995l). iGlazebrook et al.1 (Il995h 
iLe Fevre et al.l (119951), ISmail et al.1 (1 19951). iMamonI 
19981). Fostmaneta.l1 (Il998l). lArnouts et all (Il999l) 



lArnouts et al.1 (l200ll) . IMetcalfe et al.l (I2001I). lYasuda et al 
( 20011) . ICapak et all (l2004 l: if band. iMobasher et al 



19861). [Gardner et all (Il993l). ISoifer et al l 
iDiorgovski et al. (|l995l ). IGlazebrook et al.l 



iMcLeod etlo] (Il995h. iGardner et all (Il996l), iHuang et all 



il997h. IMoustakas et al.l (1 19971). ISaracco et al.1 lll997fK 
iBer shadv et al.1 (119981)7" ISzokolv et al.1 (Il998l). iMamonI 



: 1998^■ iMcCracken et al] (l2000al ). IVaisanen et al.l hOOtj], 



Martini ('2001'). 'Hu ang et ai] (l200ll). iKochanek et al .1 



(2001), Kiimmel fc Wag nerl (I2OO1I), iMaihara et al.1 (|200ll ), 
Metcalfe ot al. (2001), S aracco et al.l (I2001I). Vandame et a l. 
(2001, preprint, astro-ph/0102300 ). Mino wa et al.1 (|2005h . 
and the 2MASS Ks counts from the Durham galaxy counts 
web pages. The majority of the above number counts were 
retrieved from the Durham galaxy counts web pagefl 



A2 Amplitude of the angular correlation function 



Measurements of 
were t aken from: 
band. 



the 
U 



Roche et al 



Koo fc Szalavl (| 19841 ), 
R band 



e an gular correlation function 
band, lEfstathiou et al] (Il99ll): B 



Roche et al 



Hi 



Efstathiou et al 



^Efstathiou et al 

■ iBrainerd et al.i (jl995i). iHud on 



1991), 



(1991), 



amerg et al.i 1 

19961) ■ I Villumsen et al.l (Il9'97l), I Woods fc Fahlman 
19971 ). iRoche fc Ealed (Il999l) IConnollv et al.l (12002 



I band. 



ICabanac et al 



Efstathiou et al.l 
ll2000l). 



( I991I). jPostman et al.l 



WilsonI (l2003l). ICoil et al 



200J); K band, iBaugh et al.l (Il996l). ICarlberg et al . 
.1997I). iRoche et al.l (|l999l ). iKiimmel fc Wagnerl (|2000l) , 
iMaller et all (|2005l ). 
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